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Abstract 


Most  uudenuatcr  trackuis^  applications  require  the  estimation  ot  tar>^et  an^le,  uiiether  in 
iizunutli  or  elevation.  In  tins  note,  the  ettect  of  intesp'atwn  time  on  the  performance  of  a 
simple  tiire:el  an^le  estimator  has  been  anah/sed  via  computer  simulation.  I'lie  simulated 
estimator  mean  is  shown  to  match  theoretical  calculations  of  the  mean  of  the  estimator 
for  time-bandwidth  products  from  5  to  lOU.  Computed  variance  results  are  also 
presented.  Some  practical  considerations  for  the  simulation  of  narrow  bandwidth 
processes  are  discussed. 
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Simulation  Analysis  of  a  Target 
Angle  Estimator 

1.  Introduction 

The  area  of  torpedo  countermeasures,  for  surface  ship  and  submarine  defence,  is  a 
current  activity  at  MRL,  A  necessary  part  of  this  study  is  the  understanding  ot  the 
effect  of  different  signal  prc>cessing  structures  on  torpedo  performance.  In  this 
technical  note  a  processing  structure  that  could  be  implemented  using  analog 
techniques  is  analysed  using  computer  simulation.  Only  target  angle  estimation  is 
considered. 

In  order  to  home  onto  targets  the  torpedoes  must  estimate  the  target  azimuth  and 
elevation  angles,  relative  to  the  torpedo  heading.  The  conventional  technique  is  to 
estimate  the  time  delay,  or  phase  shift,  between  signals  received  from  two  spatially 
separated  receivers.  These  signals  would  usually  be  the  left  and  right  beam  inputs 
formed  by  summing  different  sets  of  elements  in  a  transducer  array  (if  considering 
azimuth  angle  estimation).  When  in  active  mode  the  signals  of  interest  are  often  in  a 
narrow  band  around  the  transmitted  frequency,  allowing  narrow  band  signal  theory 
to  be  used.  As  most  torpedoes  use  the  same  input  processing  for  both  active  and 
passive  mode  operation  the  narrow  band  assumption  is  also  valid  for  passive  mode 
and  is  the  approach  taken  in  this  note. 


2.  Estimator  Structure 

For  convenience  consider  the  problem  of  a  target  manoeuvring  in  a  horizontal  plane. 
The  geometry  of  the  problem  is  illustrated  in  Figure  1.  Target  physical  angle  d  is 
measured  relative  to  the  centre  of  the  array,  wlule  left  and  right  indicate  the  effective 
spatial  location  of  the  left  and  right  beam  signals  formed  by  summing  different 
elements  in  the  transducer  array.  The  target  signal  and  the  background  noise  are 
assumed  to  be  independent  white  Gaussian  noise  processes.  The  electrical  phase 


difference  between  the  left  and  right  signals  is  6,  and  is  related  to  the  target  physical 
angle  ())  by  the  transducer  beam  patterns. 


Target 


Arrav  Plane 


Figure  1:  T raiisdtiicr  array  geometry 


A  typical  input  processing  scheme  is  shown  in  Figure  2.  The  signals  are  filtered  be 
bandpass  filters,  usually  centred  at  the  receiving  frequency  for  a  stationarv  target 
(assuming  an  active  sonar),  with  the  bandwidth  covering  the  expected  range  of 
Doppler  shifts.  For  example,  a  system  with  a  transmission  frequency  of  25  kHz  ould 
require  a  2  kHz  bandwidth  to  detect  targets  with  speeds  of  up  to  ±  58  knots,  hence  the 
noise  processes  forming  the  left  and  right  input  signals,  after  the  bandpass  niters,  can 
be  considered  to  be  narrow  hand  processes.  The  scheme  as  shown  v\  ould  be  suitable 
for  passive  systems,  active  systems  could  use  the  scheme  in  con|unction  with 
appropriate  signal  detection  processing  to  indicate  whether  a  valid  return  is  present 


Figure  2:  Input  processing  scheme 


The  sum  and  difference  signals  are  formed  to  aid  the  recovery  of  the  phase  shift 
between  the  left  and  right  signals.  As  shown  in  Figure  2,  further  phase  shifting  and 
combining  of  the  sum  and  difference  signals  is  performed  before  the  signals  arc 
magnitude  detected  and  integrated  with  a  smoothing  time  constant  T,  The  minimum 


6 


integration  time  is  the  inverse  of  the  bandpass  filter  bandwidth  B  (Hz),  giving  a 
minimum  time-bandwidth  product  (BT)  of  unity.  An  estimator  is  proposed  which 
uses  the  integrated  signals  to  recover  the  electrical  phase  angle  via  the  transformation 


where 


:  2arctan(c) 


E(:b  +  fli)-  E(lf)-ai) 

In  the  above  equation  c  is  the  estimator  for  the  electrical  phase  angle  0  and  E(.v) 
stands  for  the  expected  value  of  If  the  estimator  is  unbiased  then  the  expected  value 
(mean)  of  the  estimator,  t,  will  be  0, 


3.  Analytic  Description 

3.1  Atialiftic  Signals 


For  the  svstem  in  Figure  2,  witli  white  ti.iussian  noise  input  processes  (passne 
operation),  the  left  and  right  signals  after  the  bandpass  filters  can  be  written  as 

i.  =  .Vii  cosco,,/  -  \  ;  MO  (0  :  -  I  cos(0)„f  0/ 2)-^  s.^  sin(  d),,,'  *  fi  2)  (aai 

and 

R  zz  cos  to,,  1  +  .V sin  io„f  -0,,  cos(WQf-e/2)  +  5,;  sintuiQl  -0/2)  ,  (3b) 

where  /V, ,  ,  Xh; ,  .V;.  and  .V ,,  are  independent  Gaussian  noise  sources,  representing 
the  background  noise,  each  with  zero  mean  and  variance  Ov*.  Similarlv, 

S,,  and  S,;  are  also  independent  white  Gaussian  noise  sources,  representing  the 
target  source  signal,  witii  zero  mean  and  variance  Ov'  while  ta,  is  the  centre  trequenev 
of  tile  bandpass  filters. 

,As  w'l  itter.  in  eqn.  (3),  the  input  equations  contain  both  em’eiope  and  earner 
components.  The  carrier  can  be  remiwed  from  the  equations  as  the  information 
content  of  the  signal  is  only  in  the  envelope  of  the  signal.  To  do  this  the  comple.x,  or 
analytic,  signal  representation  is  used.  If  .T(f)  is  a  real  narrow  band  signal,  its  anak  tic 
equivalent  is  denoted  .v(/)  and  is  related  to  the  real  signal  by  (IJ 


v(f)=2(x(f)e-'“'''),,. 

(4a) 

l.PF 

Y(t)=  Re!i-(f)c'“»' ]  , 

(4b) 

where  LPF  stands  for  Low  Pass  Filter,  and  j  =  v -1 

Using  this  notation,  the  input  signal  equations  can  be  rewritten  as 


and 


L  =  Se'^^+N, 

i3aJ 

r=Sc-'‘'^2  +N2  , 

(pb) 

where  S  =  S|i  =•'^'11  ^2  -‘'^'21  ~J‘^22 

In  the  next  section  tins  formulation  is  used  to  develop  the  estimator  equations  in 
terms  of  the  input  noise  processes. 


3.2  Estimator  Equations 

The  processing  structure  oi  Figure  2  is  still  valid  tor  tne  analytic  signal  represeiu.itu'n, 
noting  that  all  processing  is  now  using  complex  signals  rather  than  real  .signais  i  he 
signals  after  phase  shifting,  can  be  written  in  terms  ot  the  L  and  /<  signais  to 

form 


and 


ri=(L  +  K)c''''’ 


(naj 


b=(L-R]e-‘’'^ 


inb) 


Using  the  expressions  in  eqn.  (6),  substituting  for  L  and  R  from  (5),  and  with  some 
manipulation,  expressions  for  the  input  signals  to  the  magnitude  detection  and 
averaging  blocks  in  the  processing  structure  can  be  obtained.  I  hese  are 


/  -  e  - 

I'i  =  2Scos  — +  N]  +  iV, 


and 


f/0  0 

b  +  H  =  -v/Zt S  cos  — +  sin  —  1(1  +  /  )  +  .N/.  +  /N;  ^ 


\z(  0 

01 

\  S  cos— - 

-sin  — 

1  1  2 

2) 

)  +  N;  +;N,  [ 


IS) 


(9) 


Before  these  expressions  can  be  used  in  the  final  estimator  equation  (eqn.  (2))  the 
magnitude  of  the  complex  number  must  be  taken  and  then  averaged  over  several 
consecutive  magnitude  values.  In  general  this  calculation  would  require  the 
probability  distributions  of  the  above  quantities,  and  is  a  multivariate  problem  when 
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more  than  a  single  magnitude  measurement  is  considered.  However,  an  asvmptotic 
expression  for  the  expected  value  of  the  estimator  mean  can  be  calculated.  The 
estimator  mean  will  approach  the  asymptotic  result  for  either  long  integration 
(averaging)  length  or  large  signal-to-noise  ratios.  Let  i  =  .r.  -  /.t,,  where  x-.  and  .Xi  are 
independent  zero  mean  white  Gaussian  noise  processes  with  equal  variance  o.v^  then 
[2,  p.l95J 

E(HI)  =  e[  J.v?  +.x:;  .  (10) 

The  expected  values  can  now  be  calculated  using  eqn.  (10)  remembering  that 
5,  N,  ,  andiVj  are  independent  complex  white  Gaussian  noise  processes  with 
twice  the  variance  of  their  components.  After  some  calculation  the  expected  value  for 
c  (eqn.  (2))  can  be  shown  to  be 


,  .  ■  it  H 
j  cos  sin 


E(c)  =  • 


’  «  B  r 

■hOv  cos  -->in  j 

\  '  2  2 


.so; 


2j8o:cos'  -4oi 
\  2  • 


(ill 


As  a  result  of  the  additive  noise,  eqn.  (11)  will  be  a  biased  estimate  of  0  for  all  input 
cases  except  the  infinite  signal-to-noise  ratio.  In  this  case  eqn.  (11)  can  be  shown  to 
reduce  to 


E(c)  =  tan^  .  (12) 

Substituting  this  a-,  eqn.  (1)  yields 

E{£)  =  e,  (1.3) 

which  is  the  expected  result  for  the  estimator.  Letting  S  =  2o,v-  and  N  =  4ov-,  eqn. 
(11)  and  eqn.  (1)  can  be  combined  to  specify  the  asymptotic  performance  of  the 
proposed  estimator  for  any  signal-to-noise  ratio.  Denoting  the  asvmptotic  re.suit  bv  C 
the  result  is 


L  =  2arctant 


H 

<  0  .  0) 

cos  —  +  sin  — 

(2  2  J 

V  N 

(  0.0)' 
cos  —  sin  - 
(2  2) 

,  4S  .0 

1  +  —  cos  - 
N  2 

(14) 
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4.  Simulation  Model 


Tlie  estimator,  as  defined  in  section  2,  was  simulated  using  the  Signal  Processing 
Workstation  software  package,  running  on  a  Sun  Sparcstation  SLC  hardware 
platform.  The  Signal  Processing  Workstation,  or  SPW,  is  a  graphical  based  snnuiatn'n 
tool  which  allows  the  user  to  quickly  develop  simulation  models  for  signal  p.ocessing 
systems.  Models  are  developed  using  basic  building  blocks,  with  the  input  protessing 
shown  schematically  in  Figure  2  translating  directly  into  the  basic  building  blocks 
The  model  is  completed  with  the  addition  of  signal  source  and  sink  blocks,  allowing 
specification  of  input  signals  and  viewing  of  output  signals. 

Analytic  signal  representation  allows  significant  computational  savings  in  the 
simulation  model.  For  a  system  operating  at  25  kHz,  with  a  2  kHz  bandwidth,  the 
minimum  sampling  rate,  using  real  signals  without  bandpass  sampling,  would  be 
52  kHz.  As  the  processing  scheme  being  considered  has  phase  shifts  between  ihe 
input  signals,  the  sampling  rate  would  have  to  be  high  enough  to  allow  the  pha^e 
shifts  to  be  implemented  as  a  time  deiav  between  the  signal  components  ot  tile  input 
signals.  If  20  samples/ cvcle  at  25  kHz  are  desired  then  a  sampling  rate  ot  spo  kl  1/  w 
required.  The  z  rt/-\  phase  shifts  in  the  scheme  also  present  difficulties  it  real  signals 
are  used.  Anoilter  concern  is  the  need  to  implement  narrow  band  white  noise  sources 
for  the  input  signals.  While  this  is  straightforward,  the  calculations  add  signiticantiv 
to  the  computational  load  and  complexity  of  the  final  model 

Basing  the  mode!  on  analytic  signals  overcomes  many  of  these  concerns  T  he  ma)or 
advantage  is  that  the  required  sampling  rate  is  dictated  by  the  envelope  bandw'idth 
and  not  by  the  carrier  frequency.  As  the  signals  arc  now  complex  the  sampling  rate 
can  be  set  equal  to  the  bandwidth,  while  phase  shifts  are  implemented  bv  a 
multiplication  of  the  signal  by  c  “  where  a  is  the  phase  shift.  For  the  scheme 
considered  in  this  note,  one  result  is  that  all  references  to  frequenev  are  removed  from 
the  model  with  each  sample  representing  a  time-bandwidth  product  (BT)  of  unitv. 

The  required  noise  sources  are  now  simply  white  noise  sources  with  no  further 
processing  required.  The  computational  effort  for  the  analytic  model  is  approximateh 
three  orders  of  magnitude  less  than  for  the  real  simulation  model. 


5.  Results 

Monte  Carlo  simulations  were  performed  for  time-bandwidth  products  of  5,  10,  25,  50 
and  100.  For  each  case  100  passes  were  used  to  evaluate  the  mean  and  variance  of  the 
estimator.  Figures  3, 4  and  5  show  the  estimator  mean  results  for  BT  =  100,  25  and  5, 
with  the  asymptotic  results  calculated  using  eqn.  (14)  included  for  comparison 
The  signal-to-noise  ratios  shown  are  10  logio  S/N.  The  BT  =  100,  and  BT  =  25  are 
very  close  to  the  asymptotic  results,  while  the  estimator  mean  for  BT  =  5  is  starting  to 
deviate  from  the  asymptotic  results,  the  degradation  being  due  to  the  short 
integration  time. 
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The  effect  of  varying  the  integration  time  on  estimator  performance  is  more  clearly 
shown  in  the  variance  results.  In  figures  6,  7  and  8  the  estimator  variance, 

E((^  - 1)^  )has  been  plotted  against  different  time-bandwidth  products  for  fixed 
signal-to-noise  ratio.  The  degradation  in  estimator  performance  is  shown  by  the  rapid 
increase  in  estimator  variance  as  the  integration  time  is  decreased. 

Figures  9,  10  and  11  show  the  output  from  Monte  Carlo  runs  for  BT  =  100,  25  and 
5,  with  0  dB  signal-to-noise  ratio.  These  show  the  estimated  angles,  as  a  function  of 
time  (sample  number),  for  a  hypothetical  target  that  has  an  initial  target  electricai 
angle  of  0°  and  increases  in  20“  steps  after  every  100  samples  for  a  total  of  1000 
samples.  The  final  target  electrical  angle  is  180°.  As  expected  the  Monte  Carlo  results 
for  BT  =  too  are  localised  around  the  expected  values,  v/hile  the  variation  m  the 
expected  values  increases  as  the  BT  product  decreases. 

The  variance  results  presented  in  this  note  are  for  the  analytic  system.  If  the 
processing  were  t('  he  implemented  as  a  real  signal  svstem  then  these  i-esults  could 
only  be  used  as  an  indication  of  the  variation  m  the  real  signal  variances  as  the 
integration  time  and  signal-to-noise  ratios  are  varied.  Quantitative  values  for  the  real 
signal  variance  results  could  be  obtained  via  a  full  real  signal  simulation,  with 
associated  computational  expense.  .Alternatively  the  system  can  be  analysed  using  a 
probabilistic  approach  (2)  to  determine  the  relationship  between  the  real  and  analytic 
signal  variance  results. 
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Figure  3:  Estimator  mean  for  BT  -  100  (lines  -  asymptotic  values,  points  -  simulation 
values). 
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Figure  4:  Estwuitor  mean  for  BT  =  25  (Hues  -  iisuiuptotic  voiues,  points  -  'iniulath'ii 
values). 
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Figure  6:  Eitnnator  ;\iriiiiuv.  ratio  ~  r>  liB. 
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Figure  7:  Estimator  variance,  signal-to-noise  ratio  -  0  dB. 
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Estimate  value  ^  Estimator  variance,  degrees  squared 
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Figure  9:  Monte  Carlo  ran/ output,  BT  =  100  (signal-to-iioise  =  0  dB). 
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Sample  number 


Figure  10:  Monte  Carlo  raw  output,  BT  =  25  (signal-to-noibe  =  0  dB). 
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Figure  11:  Monte  carlo  raw  output,  BT  -  5  (signal-to-noise  -  0  dB). 


6.  Conclusion 


The  effect  of  integration  time  on  the  performance  of  a  target  angle  estimator  has  been 
analysed  via  computer  simulation.  An  asymptotic  expression  was  developed  tor  the 
estimator  mean  and  compared  to  the  simulation  results.  The  results  obtained  tor  the 
mean  are  valid  for  both  real  and  analytic  models  of  the  estimator  processing 
structure,  while  the  variance  results  can  only  be  used  as  a  qualitatii  >•  measure  ot  the 
variance  for  the  real  simulation  model. 
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